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We study the field profile of a scalar field (fi that couples to a matter fluid (dubbed a chameleon 
field) in the relativistic gravitational background of a spherically symmetric spacetime. Employing 
a linear expansion in terms of the gravitational potential $c at the surface of a compact object with 
a constant density, we derive the thin-shell field profile both inside and outside the object, as well as 
the resulting effective coupling with matter, analytically. We also carry out numerical simulations for 
the class of inverse power-law potentials V{4>) — M*+"(/)~" by employing the information provided 
I by our analytical solutions to set the boundary conditions around the centre of the object and 

show that thin-shell solutions in fact exist if the gravitational potential "I>c is smaller than 0.3, 
which marginally covers the case of neutron stars. Thus the chameleon mechanism is present in the 
relativistic gravitational backgrounds, capable of reducing the effective coupling. Since thin-shell 
, solutions are sensitive to the choice of boundary conditions, our analytic field profile is very helpful 

^■f^ ■ to provide appropriate boundary conditions for $c < 0(0. 1). 
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^ • I. INTRODUCTION 

The origin of the so called dark energy responsible for the present cosmic acceleration remains a great mystery. 
CN| ■ Since the cosmological constant (originating from the vacuum energy) is plagued by a severe fine-tuning problem, 
^ ' many alternative models have been proposed to account for the origin of dark energy (sec Refs. [l| for reviews). A 
. number of these models, including the quintessence k-essence Q and tachyon 0| models, make use of a scalar 
■ field with a very light mass (m^ ~ 10""^^ eV) in order to account for the present cosmic acceleration. If the scalar 
, field originates from candidate theories for fundamental interactions such as string theory or supergravity, it should 
interact with the standard model particles with a long ranged force (the so called "fifth force" ) . In string theory, for 
example, a dilaton field universally couples to matter as well as gravity Q . Similarly, in modified gravity theories such 
as f{R) gravity ^ and scalar-tensor theories Q, the scalar degree of freedom interacts with the matter fluid (except 
I for radiation). This is clearly seen if one transforms the action to the Einstein frame via a conformal transformation 
Q. For example, it is known that Brans-Dickc theory Q (a historically irnportant class of scalar-tensor theories) 
gives rise to a constant coupling Q between the scalar field and the matter [10]. In this sense such modified gravity 
, theories can be regarded as a coupled quintessence scenario |ll| in the Einstein frame. 
' In the absence of a scalar-field potential, the present solar-system tests constrain the strength of the coupling Q 
. 5^ I to be smaller than the order of 10"'^ [l^l- However, the couplings that appear in string theory Q and f{R) gravity 
[T^ l are typically of the order of unity. In such cases it is not possible to satisfy the local gravity constraints, unless a 
scalar-field potential with a large mass exists to suppress the coupling in the regions of high density. Moreover, if the 
same field is responsible for the cosmic acceleration today, the potential needs to be sufficiently fiat in the regions of 
low density (i.e., on cosmological scales). 

In spite of the above requirements it is possible for the large coupling models to satisfy the local gravity constraints 
through the chameleon mechanism [3, H, while at the same time for the field to have sufficiently small mass 
to lead to the present cosmic acceleration. The existence of a matter coupling gives rise to an extremum of the 
scalar-field potential around which the field can be stabilized. In high density regions, such as the interiors of the 
astrophysical objects, the field mass about the extremum would be sufficiently large to avoid the propagation of 
the fifth force. Meanwhile, the field would have a much lighter mass in the low-density environments, far away 



'Electronic address: 
t Electronic address: 
t Electronic address: 



'shinji@rs.kagu.t us.ac.jp] 

tamaki@ gravity.phys.waseda.ac.jp] 
]r tavakol@qmul.ac.uk| 



2 



from compact objects, so that it could be responsible for the present cosmic acceleration. In the case of inverse 
power-law potentials V{(j)) = M^+"(/)~" [l^ with n > 1, local gravity constraints can be satisfied for M < 10~^ eV 
[ij]- Interestingly, this roughly corresponds to the energy scale required for the cosmic acceleration today. See 
Refs. [H [O, H lli, ill, El, m, for works concerning a number of interesting aspects of the chameleon 
mechanism. 

So far the analyses of the chameleon mechanism have typically concentrated on the weak gravity backgrounds where 
the spherically symmetric metric is described by a Minkowski s pac etimc. This amounts to neglecting the backreaction 
of gravitational potential on the scalar-field equation. In Ref. [2J| the field profile in the Minkowski background was 
analytically derived both inside and outside the object by taking into account the mass of the chameleon field inside 
the body. In this settings it has been shown that the field would need to be extremely close to the maximum of the 
effective potential around the centre of the spherically symmetric body in order to allow thin-shell solutions required 
for consistency with the local gravity constraints. 

If we take into account the backreaction of gravitational potential to the field equation, the relativistic pressure 
is present even in weak gravity backgrounds such as the Sun or the Earth. It is expected that this effect changes 
the field profile inside the body in order to allow the existence of thin-shell solutions. We shall analytically derive 
the thin-shcU field profile using a linear expansion in terms of the gravitational potential $c 1) at the surface of 
compact objects. In fact wc show that there exists a region around the centre of the massive objects in which the field 
evolves toward the maximum of the effective potential because of the presence of the relativistic pressure. In order 
to realize thin-shell solutions, the driving force along the potential needs to dominate over the pressure for distances 
larger than a critical value r = r^. This distance (r^) is required to be smaller than the distance ri at which the field 
enters a thin-shell regime. In spite of such different properties of the field profile inside the body relative to the case of 
the Minkowski background, the effective coupling Qcs outside the body can be reduced by the presence of thin-shell 
solutions. We confirm this by using numerical simulations for a class of potentials of the form V{4>) = M'^~^"(j)~'^ . 

To study the viability of theories with large couplings, it is important to determine whether thin-shell solutions 
can also exist in strong gravitational backgrounds with $c ^ 0(0.1). We shall derive analytic solutions using linear 
expansions in terms of $c and then carry out numerical simulations to confirm the validity of solutions in the regimes 
with <I>c < 0(0. 1). Our analytic solutions arc useful as a way of finding the boundary conditions around the centre 
of the object in order to obtain thin-shell solutions. By choosing boundary conditions with field values larger than 
those estimated by the analytic solutions, we shall demonstrate numerically that the thin-shell solutions are present 
for backgrounds with gravitational potentials satisfying $c ^ 0.3, in the case of the field potentials of the type 
V{(j)) = A/ This marginally covers the case of neutron stars. In backgrounds with still larger gravitational 
potentials the relativistic pressure around the centre of the object is so strong that the field typically overshoots the 
maximum of the effective potential to reach the singularity at </) = 0, unless the boundary conditions of the field around 
the centre of the body are chosen to be far from the maximum of the effective potential. This overshoot behaviour is 
similar to the one recently found by Kobayashi and Maeda [2^ in the context of f{R) dark energy models (see also 
Ref. nil). We note, however, that our analytic solutions based on the linear expansion of <I>c do not cover the field 
profiles for the really strong gravitational backgrounds with $c = 0(1). In such cases we need a separate analysis 
which incorporates the formation of black holes. 

The outline of the paper is as follows. In Section [TT] wc discuss our theoretical set up as well as giving the relevant 
equations for the case of a spherically symmetric central body. In section [TTTl we give the analytical thin-shell solutions 
to the scalar field equations, both inside and outside of the body, and consider in turn the matching of thin-shell 
solutions. In Section [TVl we study the analytical field profile in more details and discuss how the field evolves as a 
function of r in the presence of the relativistic pressure. In Section |V] we integrate the field equation numerically and 
show the existence of thin-shell solutions for $c ^ 0.3. Finally Section IVll contains our conclusions. 



II. SETUP 



We consider settings in which a scalar field (f> with potential V^cj)) couples to a matter with a Lagrangian density 
Cm- In particular we shall study theories based on the action 



S = J d-^xy^ 



d4x£„(vI/W,gW) 



(1) 



where g is the determinant of the metric 5^^, A/pi = I/^/StKJ is the reduced Planck mass (G is the gravitational 
constant), i? is a Ricci scalar, and are matter fields that couple to a metric gl^l related with the Einstein frame 
metric g^i/ via 

(2) 
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Here Qi are the strength of couphngs for each matter field. In the fonowing we shall consider case s in w hich the 
couplings are the same for each matter component, i.e., Qi = Q, and use units such that Afpj = I/VSttG = 1. We 
restore G when it is needed. 

An example of a scalar-tensor theory which gives rise to constant couplings Q in the Einstein frame is given by the 
action 



S = / d^Tv/-5 



(3) 



where a tilde represents quantities in the Jordan frame. The action ([3|) is equivalent to that in Brans-Dicke theory 
with a potential C/(0). Under the conformal transformation, g^i, ~ ^fJQ obtain the action ^ in the Einstein 

frame, together with the field potential Vi^cj)) = U{4>) e^^'^ . Clearly the metric g'lll in Eq. ([2]) corresponds to the metric 
(jfj^u in the Jordan frame. 

To study chameleon fields in the relativistic gravitational background of a spherically symmetric body, we consider 
the following spherically symmetric static metric in the Einstein frame: 



d,s2 = -e2*('-)dt2 



dr^ 



1 ■ 1 
r sm 



(4) 



where ^(r) and $(r) are functions of the distance r from the centre of symmetry. For the action ([T]) the energy 
momentum tensors for the scalar field (f> and the matter are given, respectively, by 



f.LV 



T. 



2 diLrn. 



-5 f^ff^" 

Under the gravitational background ([1]), the (00) and (11) components for the energy momentum tensors are 

1 



-2* i,/2 



y^-y(<^), t; 

where a prime represents a derivative with respect to r and 

rpO(m) _ g4Q0j.O(m) rpl{m) _ ^4Qtj,rpr(m) 



(5) 
(6) 

(7) 

(8) 



Here f "("^ and f^^"''^ are the energy momentum tensors of matter in the Jordan frame. Denoting the energy density 
and the pressure of the matter in the Jordan frame as /5,„ and Pm, the matter energy- momentum tensor in this frame 
takes the form T^f ~ {~Pm,Pm,PrrnPm)- The corresponding expressions for the energy density and pressure in the 
Einstein frame are then given by p,n = e^^'^pm and p„ 
The evolution equation for the scalar field (j) is given by 



dx^ d{d(f)/dx^) 



0, 



where the derivative of £„ 



^niigp.^) = ^mie^'^'gfiiy) in terms of is 



dCr, 



-gQe^^^~gf.uf^'' = V^Qi-Pm + ipm) ■ 



(9) 



(10) 



We then obtain 



1 



„2'I' 



The Einstein equations give: 



2r 

32* -1 



2r 



AnGr 



Xj/' — <j)' 



^1 



= -SttG 



^0'2 + e2*y(</))-e2*p„ 



(11) 

(12) 
(13) 
(14) 
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From the conservation equation, V^Tl* = 0, we also obtain 

V'm + (Pm + Pm)*' + (90'(Pm " 3p„) = , (15) 

which is the generahzation of the Tohnan-Oppenheimer-Volkoff equation. Note that this equation can also be derived 
by combining Eqs. pTjl - fH)) . 

Our main interest is the case in which the field potential ^(0) is responsible for dark energy. In that case both 
and 0'^ are negligible relative to pm hi the local regions whose density is much larger than the cosmological one 
(po 10-29 g/cm^). Then Eq. ^ can be integrated to give 



g2<E.(r) 



m{r) = / 47rr'2pm dr' . 



(16) 



Substituting Eqs. HI]) and into Eq. HIl) gives 

47rG're2* (p,„ - p.„ 



l + e2* 



0' = e^* [y^ + 0(p,„ - 3p„0] . (17) 



We assume that the energy density is constant inside (p^ — Pa) and outside = Ps) of the spherically symmetric 
body with a radius re- Strictly speaking the conserved density pm in the Einstein frame is given by pm = e~^'^ p^ [isl. 
QIll. However, since the condition ^ 1 holds in most cases of interest, we do not need to distinguish between 

Pm and Pm.. 

Inside the spherically symmetric body (0 < r < Tc) we have mir) = Airr^ pa/3 and Eq. (jl6[) gives 

e^*M = (l-^p^r2)"\ (18) 

With the neglect of the scalar-field contributions in Eqs. (|12p - p^ it is known that the background gravitational field 
for < r < Tc corresponds to the Schwarzschild interior solution. In this case the pressure Pmii") inside the body 
relative to the density pA can be analytically expressed as 



,(r) _ VI - 2(rVr2)$c - VI - 2$c 



(0 < r < r,) , (19) 



PA 3Vl~^2$; - Vl-2(r2/r-2)4>c 
where $c is the gravitational potential at the surface of body: 

Vq d 

Here Ale = 4:TTrlpA/i is the mass of the spherically symmetric body, and in the last equality in Eq. we have used 
units such that G = I/Stt. Equation (fTO|) shows that the pressure vanishes at the surface of the body {pm{rc) = 0). 

In the following we shall derive analytic solutions for Eq. ([TT|) . under the conditions |<I>(r)| ^ 1 and |^(?')| ^ 1- 
We neglect the terms higher than the linear order in <I>(r) and ^'(r). From Eqs. P^ - ipO)) it then follows that 

$(r) c , ^ ^ i£ A _ for < r < . (21) 

PA 2 \ riJ 

At the centre of the body we have Pm {0)/pA — which shows that the effect of the pressure becomes important 

in strong gravitational backgrounds. 

Outside the body we assume that the density ps is very much smaller than pA with a vanishing pressure. Then 
the metric outside the body can be approximated by the Schwarzschild exterior solution: 

, , GAI Tc , , , , 

$(r)~ — = p™(r)~0 for r>r,. (22) 



III. MATCHING SOLUTIONS OF THE CHAMELEON SCALAR FIELD 



In this section we solve the scalar-field equation (fTT)) in the relativistic gravitational backgrounds discussed in 

secnn 
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In the nonrelativistic gravitational background where the pressure Pm as well as the gravitational potential <I>c are 
negligible, the effective potential for the scalar field is defined as [H, IT4 | 

V,s{4>) = V{<t>)+Qp„,(l>. (23) 

This potential has a minimum either when (i) < and Q > or (ii) V_tt> > and Q < 0. An example of class 
of potentials satisfying (i) is provided by the inverse power-law potentials V{(f)) = A/'*+"0~" (n > 0). Since f{R) 
gravity corresponds to the coupling Q = — 1/\/6, the effective potential VcS has a minimum for the case > (as 
in the case of the models proposed in Refs. [23, [2^, 2Be l30i, iMi, .32.] ) . 

For constant matter densities, pA and pB, inside and outside of the body, the effective potential (|23|) has two minima 
at the field values (j)A and (j)B characterized by the conditions 

= 0, (24) 
V^i(j)B) + QpB = 0. (25) 

The former corresponds to the region with a high density (interior of the body) that gives rise to a heavy mass squared 
a), whereas the latter corresponds to the lower density region (exterior of the body) with a lighter mass 

b). 



2 _ d^Veff 



squared nig = '^^f ( 



The following boundary conditions are imposed at r = and r ^ oo: 

^(r = 0) = 0, 0(r^oo) = (/) 
dr 



(26) 



We need to consider the potential (— VJ^ff) in order to find the "dynamics" of (j) with respect to r. This means that 
the effective potential (— T4fr) has a maximum at = (f>A- The field is at rest at r = and begins to roll down the 
potential when the matter-coupling term QpA becomes important at a radius ri. If the field value at r = is close 
to (j)A , the field stays around (pA h^ the region < r < ri . The body has a thin-shell if ri is close to the radius Tc of 
the body. 

The position of the minimum given in Eq. (|24p is shifted in the relativistic gravitational background. In the following 
we shall derive the field profile by taking into account the corrections coming from the gravitational potential. Inside 
the body, Eq. ([T7]) to the the linear order in $c reduces to: 



1 $ 

2rr' 



0' - (1<0 + Qpa) 1 + 2$, 



QpA<^c 1 - — = 



(27) 



In the region < ?' < ri the field derivative of the effective potential around (j)A may be approximated by 
dVcff/d4> = + Qpa — ^1^(0 — (I)a)- The solution to Eq. (P?)) can be obtained by writing the field as 4) = (pQ + S(j), 



where (pQ is the solution in the Minkowski background and 
in S(j) and $c we obtain 



is the perturbation induced by At the linear order 



+ -00 - "i^('/'o - 0a) = , 



^QPA 



The solution to Eq. that is regular at r = is given by (poir) = (j)A + A{e 
Substituting this solution into Eq. ((29|) we obtain the following solution for p{r): 



(28) 
(29) 



e"^^^)/r, where A is a constant. 



(r) = (f>A + 



A{e- 



1 



rriAr' 



3QpA<^c 



1 2 2 1 



1 



-TJiAr - - 
4 4 



8m Ar J \3 
(0 < r < 7'i) . 



1 2 2 . 1 
—rriAr 



—mAr 

4 4 



1 



8m^r 



(30) 



One can easily show that this solution satisfies the first of the boundary conditions ([26|) . 

In the region ri < r < the field |0(?')| evolves towards larger values with increasing r. Since jV^^j <^ \Qpa\ in 
this regime one has dVcs/dcl) ~ Qpa- In this case (pQ and Sep satisfy 

00 + l<f>'o -QPA^O, (31) 

„2 - 



-Sep' = $c 
r 



^/o-IQpaU-I^- 



(32) 
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Wc then find the foUowing solution 



B f \ 1 , / 3 23 



(ri < r < Tc) 



(33) 



where S and C are constants. 

The field acquires suflicient kinetic energy in the thin-shell regime, in order to allow it to climb up the potential 
hill towards larger absolute values in the region outside the body. As long as the kinetic energy of the field dominates 
over its potential energy, the right hand side of Eq. (|17p can be neglected relative to its left hand side. Also the term 
that includes pm and pm in the square bracket on the left hand side of Eq. ()17p can be neglected relative to the term 



(1 



o2* 



)/r. Using Eq. 



the field equation reduces to 



The solution to this equation is 



(r) = + 



D 



GM 



GM 



(r > rc) 



(34) 



(35) 



where Z? is a constant. Note that here we have used the second boundary condition in Eq. (pS)) . 

Having obtained the solutions ([50]) . ([55]) and (|35p in the three regions inside and outside the central body, we 
proceed to match these solutions at r = ri and r ^ Vc- The thin-shell corresponds to the region defined by 



Arc = rc - ri < Tc 



(36) 



namely Arc/rc ^ 1. 

It is possible to satisfy the local gravity constraints as long as the field inside the body is sufficiently massive, i.e., 
"m-ATc ^ 1 (or niAri ^ 1). For example, in the case of the Earth with the class of potentials V{(f>) = M^+"</)~", 
we have the constraints mATc ^ 10^, for n = 1, and mATc ^ 10^, for n = 2, from the experimental tests of the 
equivalence principle [2^ . Since the mass itia becomes larger in higher density regions, the quantity m^T'c inside a 
strong gravitational body becomes even larger than in the case of the Sun or the Earth. We use the approximation 



that 



is negligible relative to e™'*''^ in Eq. ([3^ 



Wc recah that $c = O{10-'^)-O{10-^) for neutron stars, $c = O{10-^)-O{10-^) for white dwarfs, $c = O(10~^) 
for the Sun and <i>c = 0(10"^) for the Earth. In the following we shall use linear expansions in terms of the three 
parameters Arc/rc, <&c and l/(?n^rc) (or l/irnAfi))- We drop terms of higher order in these parameters relative 
to 1. We caution that our approximation loses its accuracy under the really strong gravitational backgrounds with 
*c>O(0.1). 

Using the continuity of 0(r) and </''(?') at r = ri and r ~ rc we obtain 



ri 
ruAe^' 



ruArl^c 



ri 



-irl 
mAr\^ 
3r2 



1 



AmATi 



1 / r\ 



ri 



+ C = 4>, 



QpArl ( 3 23 ?i \ 



AmAT\ 



TfiArx 



■2rl 



1 + ^cTT^ = -T-QPATi 2 - 3$c + ^*c4 . 
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23 



Bl(l-|)-C + i.l(l + <I>.)^ 



B\\ 



+ /?(! + 2$,) 



(38) 
(39) 

(40) 



The value of C can be derived from Eqs. ([55]) and keeping terms to linear order in (f>c only. Substituting C 
into Eq. ((37)) and using Eq. p8|) wc can obtain expressions for A and B. The coefficient D is then obtained from 
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Eq. ([40]) . Using this procedure we find 



A 



1 



QpArf 1 - -$c 



3r2 



7^ ^? 

4 ^2 



1-$, 



'2r2 



A~cf>B){l + <^c^]+:^QpArUl- 



'=2r2; 
(41) 



B ^ -{l-a)ri 



9. 



7^ 

4^r2 



1 



D = (1 - a)ri 
1 



Qp^r^ ( 1 - + $,2^ ) - -QpArf ( 1 - 3$, 



6 



1 



QpAvi 1 - J + -Qp^r-f 1-3$ 



5 V2 



where 



(r2/3r2)$, + l/(mAri) 



1 + (r2/4r2)$, + (?7iArf$c/3r2)(l - {rf/2r^^)^c) 



Since the denominator in Eq. (j45p is larger than 1. the parameter a is much smaller than 1. 
The distance ri is determined by the condition 



where 



(f>{ri) = (f>A~A- 



ri 



'"A [^(ri) - 0a] = QpA : 



Substituting Eq. (gT]) into Eq. (gH) gives 

A = - 

From Eqs. (|4T|) and ([48|) we then obtain 

0A - 0s = -QpatI 
where 

/3 



3r2 



AmATi 



3r2 



4 r2 



Arc 



1 Ar,. 
2 



■mATc 



1-^1(1-/3) 



(mAr^$c/3rg)(r2/rg)$c 
1 + (mAJ-f $c/3r2) _ (r2/4r2)$c 



Note that /3 < 1. 

The thin-shell parameter introduced in Refs. [H, [l3| is in this case given by 



(42) 
(43) 



4 'r? 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



6g$ 

Ar 



1 Arc 

2 re 



m^rc 



1-^1(1-/3). 



(51) 
(52) 
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To the first-order in expansion parameters one has eth = ^fc/fc + 1/('tt-a''c)j which is identical to the corresponding 
value derived in the Minkowski background [IJ]- The effect of the gravitational potential appears as a second-order 
term to the thin-shell parameter. Substituting Eq. (j49ll into Eq. ((44|) . we obtain the following approximate solution 



D ~ -QQ^cTc 



An 



1 



Arc 



1 



rriATc 



Arr 

1-2 --<S>c-a-fi 



(53) 



where we have carried out a linear expansion in terms of a, /3, Ar^jr^ and $c- The solution outside the body is then 
given by 



(r) ~ 0s - 2Qeff 



GM 



1 



GM\ 



where the effective coupling is 

Qeff = 3(5 



Are 



Arc 
rr 



m^rc 



Ar 



(54) 



(55) 



To leading-order this gives Qeff = 3Q(Arc/rc -I- l/m^rc) = 3(5eth, which agrees with the corresponding result in the 
Minkowski background [l^]. Thus provided that eth ^ Ij the effective coupling Qgff becomes much smaller than the 
bare coupling Q. The gravitational potential <I>c appears as a next-order term. As can be seen from Eq. (j55p the 
presence of the gravitational potential <i>c leads to a small decrease in Qofr compared to the nonrelativistic gravitational 
background. 



IV. THE FIELD PROFILE 



In this section we shall discuss the analytical field profile derived in the previous section in more details. The 
coefficient A and the field difference c^a — <i>B are determined by fixing the value of ri, see Eqs. and From 
Eqs. dSO]), (iMl), dSSl), gU-dlil), dm) and jMl) the thin-shell field profile is given by 



0(r) 
where 

B E 



^cri QPA 



mArl<Pc <^cri 



,2 „2 



-rriAr 



1 1 

4 SjTi^r 



4^2 

^' c 

4r2 



(e" 



1 + 



2m\ 



6 



0(r) =0^-1- 



6eth + 6B— 
r 



1 - 



2^ 



r V r 



Qpat'cTi 



-eth 1 



rl f 3 



D 



QpatI 



(1-a) 



eth- 



^"1 



" 2r2 
9<I>cr^ 
5r2 

-I L_L 

2r2 



6 \ r? 



1 2 2^1 



-3 1- 



1 / *c 

+ 1 - ^ 



3$, 
~2" 



1 - 3$,. 



2re 
Sr^ 



1 



Smyir 



(0 < r < ri), 



(56) 



1 

2 



20r2 



(ri < r < rc), 



(57) 
(58) 



2r2 



2r2 



2 



2r? 



2r3 



4r? 



1 - 3$e + 



(59) 



4r2 



2\ 1 



(60) 



Note that the field profile given in Eqs. ([56)l - (|58p has been derived without specifying the form of the potential. While 
the term pA in Eqs. can be replaced by 6<i>c/r2, we have chosen not to do this so that the field profile in the 

Minkowski background can be simply recovered by setting $c = 0. 
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In the following we shall consider in details the field profile in three regions: (i) < r < r2 = 1/m/i, (ii) r2 < r < ri 
and (iii) r > ri, and discuss the case Q > for simplicity. 



A. The region < r < r2 



Deep inside the body where the distance r satisfies the condition r <C Eq. ([55]) gives the following approximate 

field value and its derivative with respect to r: 



[r) ~(I>A + 



(j)'{r) ~ QpAT, 



2QpAri 
2mAri 



1 



3r2 4r2 



6 2(mAr-c)2 



2m^ 



6 



(mArc)2 



3r? 



3*c 



r 



(61) 
(62) 



In the Minkowski background ($c — 0) we have 0(0) :^ (pA + 2QpAri/ {mAC™-^'^^) and 0'(r) > (where r 7^ 0). 
Hence the field rolls down the potential toward larger with increasing r. In the presence of the gravitational 
potential $c, the derivative (j)' {r) can be negative depending on model parameters. Using the approximation Tc — ri 
the condition that (t>'{r) < translates into 

(mAri)3 3^2 9* 
gmATi 2 ^ 2 V 4 



(63) 



When <f>c > 0.253 this is automatically satisfied for all (positive) niAri- When $c = 10^^, 10~^, 10~^, the condition 
is satisfied for niAri > 6, ruAri > 22 and niAri > 29, respectively. Hence in most realistic cases where m^ri » 1 
we have 4>'{r) < in the region < r < ?'2. Interestingly this property persists even in the weak gravitational 
backgrounds, such as those of the Sun or the Earth. 

The evolution towards the smaller (f> region is due to the effects of the relativistic pressure Pm since the last term 
on the right hand side of Eq. (pT|) originates from the pressure. Compared to the case of the Minkowski spacetime, 
the presence of the pressure term leads to the shift of the field (j){Q) towards a larger value. When this effect of the 
pressure dominates over the rolling down effect along the potential, we have 



This shows that the field decreases from (j){Q) ~ ( 
with increasing r. 



1-^ 



6 



[1 



{niArc)^ 
G/imArc)"^] 



to 0(r2) 



3Qp4*c 
2m'i 



[1 



(64) 

7/(TOArc)'] 



The region 7-2 < r < ri 



In the region r2 < r < ri we have the following approximate solutions from Eq. ([5; 
<l>'{r) 



(1>A 


+ 


Qpa mA(r-ri) 


ri/r + niA'^crir^ /ir^. 


- $crir/4r2 
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1 + ruAr^'^c/i'i'l — 




Qpa 


/ mArl<^>c 


-*;?) 


n) 


X 




<^cJnAfif^ 


^cTir 1 / ri 


2(^c'niArir'^ 


r 


3r2 


4r2 ruAr \ r 





2m\ 



6 



[mArcY 



3QpA^c 
(mA^c)^ 



(65) 



(66) 



In Eq. ((SS]) we have taken into account the terms {m\r'^ /3 — niAr /4:)e™^^ in the last square bracket of Eq. (|56p . If 
4>'{r) < at r = r2, the field derivative needs to change its sign from negative to positive at the distance (i-C, 
f (ra) = 0). 

From Eq. the condition (()'{r^) = translates to 



rriArce 



iA{ra-ri) 



1 



= 3$,^ ( 1 





T3 




4r2 iTiAr^ 


mArl<^c 








3r2 


4r2 ) ■ 





3r2 



4r2 



(67) 
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rriAri = 10 


mATi = 10"^ 


ruAri = W 


$c = 10"^ 


(?!>'(r) > 


0.75 


0.97 


$c = 10"'^ 


(f>'{r) > 


0.82 


0.98 


$c = 10"' 


0.55 


0.94 


0.99 



Table I: The values of rs/ri at which (jj'ir^) = under the approximation rc — r\. It is clear that ra/ri increases for larger 
$c and iriAri. In the cases "I>c = 10~^ and = 10"*' with niAri = 10 the field derivatives (;/>'(r) are positive in the region 
l/niA < r < ri. 



In Tabic I wc show the values of r^/ri under the approximation ~ ri for several different choices of a-nd m^ri. 
Clearly rs/ri gets larger with increasing (f>c and niAri. When (f>c = 10~^ and <i>c = 10"^ with m^ri = 10 the field 
satisfies the condition 4>'{r) > in the region r2 < r < ri. Meanwhile, when (f>c = 10"^ and toa'T = 10, the sign 
change of ^'(r) occurs at r^/ri = 0.55. 

In the region r > the second term on the right hand side of Eq. (|65p dominates over the third one, giving the 
following solution 

At r = ri we have 4>{ri) = (j)A + QpA/''n\, as required by Eq. (^5]) . Note that — (j)A\ is larger than |0(O) — (j)A\ 

by a factor of 2/(3$c). 



C. The region r > ri 

In the region j'l < r < rc the field grows because of the dominance of the last term in Eq. (j57p . The field value 
at the surface of the body can be estimated as 



0(rc) ~ + QpArl 



1 



mATc 



Arc 



1 Ar 

rUATc Tc 



m , 



1 

^+2 



ruArc- 



, (69) 



where in the second approximate equality we have used the fact that a is of the order of 1/ [niAri). Obviously 4>{rc) 
is larger than 4>{ri). If the condition 



rc ruArc 



(70) 



is satisfied, we have that (j){rc) 3> (/'(t'i). In this case the field acquires a sufficient amount of kinetic energy so that 
the following condition is satisfied dX r ~ rc'. 



QPA > \V4 



(71) 



We recall that outside the body the density pA sharply drops down to pB- Hence only the potential-dependent term 
remains on the right hand side of Eq. ([TT]). Under the condition ([7T|) the kinetic energy dominates over \V^\ for r > rc 
so that the field equation is approximately given by Eq. ([34| . In this case the analytic field profile should be trustable. 

If the condition (|70p is not satisfied, the field (j){rc) is not much different from 4>{ri). In this case the kinetic energy 
of the field is not sufficiently large so that the term \V,^\ is not negligible relative to QpA ai r = rc- In the region 
r > Tc, this can lead to the puUback of the field because the kinetic energy is not large enough for the field to climb 
up the potential hill. In fact we have numerically confirmed this behaviour in cases where Ardrc is smaller than the 
order of XjiniArc)- Thus the condition (j70p is important in order to obtain the field solution (j35p outside the body. 



V. NUMERICAL SIMULATIONS 



In this section we shall numerically confirm the analytic field profile presented in the previous section and discuss 
the validity of the approximations used to derive it. In these numerical simulations wc employ the class of inverse 
power-law potentials 

V{^) = m4+"(/)-" (n > 0) . (72) 
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Although we specify the field potential to be of the form (l72l) . the thin-shell field profile we will derive numercially 
in this section also holds for other potentials, such as V{(l)) = M** exp(Af"/0"), as can be expected from the general 
form of the Eqs. ([56|) -([58 |) . The effective potential Vee has extrema inside and outside the body for Q > 0. The field 
value (f>A and the mass squared inside the body are given by 



l/(n+l) 



_Q PA 

which lead to the following relation 



TO^ = + 1) [ Pa 



{n + l)QpArl 



(ri+2)/(n+l) 
V \ ^_(„+4)/(„+l) 



{mATcY 



The field value (pB and the mass tob can be obtained by replacing pA for pB in Eq. ([73 
We introduce a dimensionless field Lp defined by 



ip = (picpA. 

From Eqs. ([55)) - (|58p the analytic thin-shell field profile for the potential ([7^ is given by 



(73) 



(74) 



(75) 



ip{r) 



1 



+ 



1 



{n + l)e"^Ari J, 
(n -I- \)e'^'Ari ruArl 



1 2 2 1 
3™A^ - -^mAT 



ri / mArl<^c _ £c fj. 

3r2 4 rl 



3$r 



2(71+1) 



1 



6 



(niArcy 



3r? 4 r^. 



SniAT 



(fir) = 1 + 



ip{r) 



1 



n + 1 

{■mATcY 
?i + 1 



~ri / 4>cH\ 1 / 1 



23$, 



20r2 



(l + $c- 

r \ r 



(r > Tc) ■ 



Introducing a dimensionless distance normalized by rc'. 

the field equations to be solved numerically take the forms 

dV ^ 2 - 5$ea;^ + 2,^cX^Pm/Pm dy> _ {mATcf 
da; 2 



1 



x(l - 2$ca;^) dx n -I- 1 1 - 2$ca;2 

dV 2(1 - $c/x) - 3$ea:Vi3/pA d(yg _ (m^rj^ 1 
dx2 ^ 



1 - 



3pm(?') 



x- 2$c 



Ax 



n+l 1 - 2<f>c/a; V PA 



PA 
PB 1 



,n+l 



(0 < r < n), 

(76) 

(ri < r < Tc), 

(77) 
(78) 



(79) 

(80) 
(81) 



(0 < X < 1) ^ 
(•^ > 1) , 



where Pmir)/ pA is given in Eq. (fT9|) . Given the occurrence of a; in the denominator of the second term of Eq. (|80p . 
the numerical solutions cannot start from the centre of the body {x ~ 0). Instead we start the integrations from a 
radius r = ri, slightly away from the centre satisfying the condition ^ I/uia- In so doing we use the analytic 
solution ([75)1 with the field derivative 



dip 
dx 



1 



n+1 



1 



1 



3r? 



2 „2 



4 r2 
1 1 



1 



^ 2 2 

-m, aT 



^ 12 2 8771^7- 877l^r2 

5 11 1 

77l4r 7^ — - 

12 ^ 2 877iAr- 877i^r2 



3$. 



(82) 
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From Eq. ^ we have that 0b/(/'a = (pa/pb)^/^"+^^ and hence eth = (ri + l)/(TOArc)^[(pA/pB)^^^""^^^ - !]• Using 
Eq. ([52]) we obtain the following relation for the ratio pa/ Pb 

1 ii+i 



PA_ 
PB 



1 + 



n + 1 



Arc 



1 



2 Tr, 



1 



mATc 



1 



^1(1 



(83) 



Thus specifying the values of n, to^Tc, l^rdvc and $c, allows the ratio pa/ PB to be determined. Alternatively, given 
the ratio paI PB together with n and <I>c, allows the relationship between mATc and Arc/rc to be derived. We note 
that the condition Arc/rc ^ l/(myirc) needs to be satisfied for the field to have a sufiicient kinetic energy outside 
the body. 



A. Minkowski background (<3?c = 0) 



Let us first consider the Minkowski background ($c = 0). In this case the analytic field profile is given by 

1 



^p{x) = 1 + 
= 1 + 
if>{x) = 1 + 



(n + l)e™^''i rc a; 
(m^rc)^ 



^1 ^ l^^niArcX 



') (0<r<ri), 



n+l 

{mArcY 
n+l 



^-+6^^ -^)+3i^-^- 



l r? 1 



mATi / rriyi?'c J^c ^ 

2 





(84) 
(85) 

(86) 



In deriving Eq. (j86p we have used the relation + [(''i/^c)^ ^ l]/2 = ^i/(™a''c) coming from Eq. ([51 

For the parameter values n = 1, = 1 g/cm^, pB = lO^'' g/cm'^ and /^rc/vc = 0.0625 used in the numerical 
simulation of Ref. [l3|, we obtain l/(myirc) = 0.0200 and eth = 0.0793 from Eqs. (|83p and ([5^ . In this case the 
condition l^Vc/rc > l/{mATc) is satisfied so that the field acquires sufficient kinetic energy in the region ri < r < rc. 
One can also consider the case in which the difference of l^Vc/rc and l/(mArc) is larger. For example, with n = \, 
PA = 1 g/cm3, PB = 2.0 X 10"^ g/cm3, /:^rc/rc = 0.08, one has l/(mArc) = 0.0142 and eth = 0.0899. In Fig. [Tj 
we plot the thin-shell field profile for this latter case by choosing the boundary conditions for </? and = dLp/dx at 
Xi = ri/vc = 10~^, using the analytic solution ([84[) . The numerical solution (a) derived by solving Eqs. ([80[l and ([8T]) 
shows fairly good agreement with the analytic solution (b) in the region r < ri = 0.92rc. On the other hand the 
agreement is not very good in the region r > ri, with a 20 % difference at the distance r = 5rc. 
Inside the body the following relation holds 

^(-i)-1 + 4t' ?T^-4t' (87) 
n + l QpA (yJ + 

which gives \V^^{ri)\/QpA — 1/2. This shows that our analytic estimation does not hold well in the region around 
r = ri. In particular the neglect of the term relative to QpA in the region ri < r < Tc gives rise to an error 
compared to the numerical simulation including this term. We find that the field value numerically obtained in the 
region ri < r < rc is smaller than the analytic value given in Eq. (j85p . This leads to the smaller field derivative ip' (x) 
at the surface of the body (.t = 1). In the numerical solution presented in Fig. [Tjthe numerical value of ip'{x = I) is 
different from its corresponding analytic value by about 18 %. This difference is inherited by the field profile outside 
the body. 

In Fig. [T] we also plot the numerical solution (c) derived by solving the field equation with the approximation 
Vce,4, = Tn\{(l) — (Pa) for < r < ri and 14ff_0 = QpA for ri < r < Tc- We find that the solution (c) agrees well with 
the analytic solution (b). This shows that the reason for the discrepancy between the solutions (a) and (b) is due to 
the fact that the matching of two analytic solutions at r = ri overestimates the field values and their derivatives in 
the region ri < r < rc. We have also tried other model parameters and have found that this property holds generally. 

If we take boundary conditions with larger values of ip{x) or (/5'(x) than those estimated by Eq. (|84p around the 
centre of the body, it is possible to obtain a field profile outside the body that is close to the analytic estimation (|86p . 
In Fig. [2] we show the numerical solution corresponding to the same model parameters as given in Fig. [Tjbut with a 
boundary condition for the field that is larger than the one given by the analytic solution ([M[) . As can be seen in this 
case the numerical solution outside the body agrees well with the analytic solution ([55]) . Note that the field approaches 
the asymptotic value (pB = 4'b/4'A = 223.6, estimated analytically using the relation ipB = {pa/ PbY^'^"'^^^ ■ 

The above results show that the analytic solution is useful to find boundary conditions in order to determine the 
thin-shell field profile. If we choose the field value to be slightly larger than the one estimated by Eq. ([M[) around 
the centre of the body, we are able to find a numerical solution outside the body that is close to the analytic solution 
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2.0 10- 




Figure 1: The thin-shell field profile in the Minkowski background for n = 1, Q = 1, pa/pb = 5.0 x 10*, and /S.rc/rc ~ 0.08. 
This case corresponds to l/(mAJ"c) = 0.0142 and eth = 0.0899. The boundary conditions for (p and (p' at Xi — 10~^ are chosen 
by using the analytic solution (|84|l . In the left panel the black curve (a) shows the numerically integrated solution. The dotted 
curve (b) corresponds to the analytic field profile given in Eqs. (|84p - (|86p . The dashed curve (c) corresponds to the numerical 
solution that is derived by solving the field equations using the approximations VcB,rf, = m\{4> — (j)A) for < r < ri and 
Vcff,(/) ~ QpA for ri < r < Vc- While the curve (c) agrees with the curve (b) with high accuracy, the curve (a) deviates from 
the curve (b) in the region r > ri = 0.92rc. This shows that the analytic estimation that connects two solutions at r = ri 
overestimates the field value outside the body (about 20 % larger at the distance r = 5rc in this case). The right panel is 
the magnified log plot of {ip — 1) in the region < r/vc < 1.4. While the numerical solution (a) agrees well with the analytic 
solution in the region r < ri, the deviation begins to appear in the region r > ri (in the log plot the deviation appears to be 
small) . 
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Figure 2: The thin-shell field profile in the Minkowski background for the same model parameters as given in Fig. [T] but with 
a different boundary condition for (p at Xi — 10~^: ip{xi) — 1 = 1.630 x 10"^*' (which is larger than the one used in Fig. [T] 
ip{xi) - 1 = 4.889 X 10'^''). The derivative (p'{x) at x ^ xt is the same as in the case of Fig. [1] (p'{xi) = 8.074 x 10~^^. The 
black curve (a) and the dotted curve (b) show the numerically integrated solution and the analytic field profile, respectively. 
The left panel is the plot in the region < r /vc < 50, whereas the right panel is the magnified log plot of {tp — 1) in the region 
< r/vc < 5. In this case the numerical solution approaches the asymptotic field value ipB = (fiB/4'A = 223.6 in the limit 
r/rc oo. In the region outside the body the analytic solution agrees well with the numerical solution. 
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Figure 3: The thin-shell field profile for = 7.0 x 10"^°, n = 3, Q = 1, Arc/r^ = 0.085 and m^rc = 40.0. This case 
corresponds to pa/ Rb = 3.3 x 10*", cfA = 1.05 x 10"^^, cI>b = 4.48 x 10"^" and eth = 0.104. The boundary condition for the 
field at X, = lO"'^ is Lp{x^) - 1 = 2.615210 x 10"^°, which is slightly larger than the analytic value tp{xi) - 1 = 2.615179 x 10"^° 
that comes from Eq. (|76|l . The derivative (p'{xi) is chosen to be the same as the analytic value. The left panel shows ip in 
the region < r/vc < 30, while the right panel depicts {ip — 1) in the region < r/rc < 5 with log scales in the vertical and 
horizontal axes. The black curve (a) shows the numerically integrated solution, while the dotted curve (b) is the analytic field 
profile given in Eqs. (|76|) - (|78| ). The solution approaches the asymptotic value (ps = 42.705. 



B. The relativistic gravitational background ($c 7^ 0) 

We shall proceed to the case of the relativistic gravitational background. As we already explained in Sec. IIVI the 
presence of a relativistic pressure is important around the centre of the body. This relativistic pressure gives rise to 
a force against the driving force that comes from the slope of the field potential. If the condition ([55)1 is satisfied, the 
field evolves toward smaller values in the region < r < r2 = 1/toa- In this case the field derivative (^'(r) needs to 
change sign at r = ra (r2 < ^3 < ri) for the realisation of the thin-shell solution. If (p'{r) is positive in the region 
< r < ri, the field dynamics is similar to the one in the Minkowski background, discussed in the previous subsection. 

Let us now consider the case < in the region < ?' < ra. As long as m^ri 3> 1 this situation naturally 

appears even in weak gravity backgrounds such as in the case of the Earth or the Sun. In Fig. Owe present an example 
of a numerically integrated field profile corresponding to the gravitational potential of the Earth with $c = 7.0 x 10""'^° 
and n = 3, Q = 1, Avc/rc = 0.085 and m^rc = 40.0. We choose the boundary condition of (f> at Xj = 10^^ to be 
slightly larger than the analytic value derived from Eq. ([75]) . so that the numerical solution approaches the field 
value ipB = 4'b/4'a = 42.705 asymptotically. The reason for this choice is that matching two analytic solutions at 
r = ri leads to an overestimation of the field value by neglecting the term V^^ relative to the term QpA in the region 
ri < r < Tc- The resulting field profile is sensitive to a slight change of the boundary condition for ip{xi). This shows 
the importance to derive analytic solutions for finding appropriate thin-shell solutions, as we have done in previous 
sections. 

In the region < r < 7*3 ~ 0.26ri the derivative (f>'{r) is negative. It changes sign at r ~ and the field begins 
to grow in the region r > r^, for increasing r. This behaviour is confirmed in the right panel of Fig. [3) Around the 
surface of the body the field acquires sufficient kinetic energy so that it climbs up the potential hill toward (j) = (f>B ■ 
The left panel of Fig. [3] shows that the numerical solution outside the body agrees well with the analytic thin-shell 
solution given in Eq. (j78[) . Thus the chameleon mechanism is present in the relativistic background with weak gravity 

($c< 1). 

For larger gravitational potential $c, the effect of the relativistic pressure becomes stronger around the centre of the 
body. This leads to the rapid evolution of the field (p toward smaller values. We also note that the analytic thin-shell 
solution ((76|) - ((78|) begins to lose its accuracy in the stronger gravitational backgrounds with <I>c > 0.1. If we run our 
numerical code by choosing boundary conditions for (j){r) and (j>'{r) around the centre of the body determined by 
Eq. ([75)) . the solutions with $c ^ 0.1 typically keep evolving toward smaller regions by overshooting the effective 
potential maximum at </> = (j)A- However, if we choose boundary conditions for (f> which are larger than the one 
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Figure 4: The thin-shell field profile for <E>c = 0.2, n — 2, Q = 1, Avc/vc — 0.1 and ruArc ~ 20.0. This case corresponds to 
Pa/pb = 1.04 X lO"*, cl)A = 8.99 x 10"^, (f)B = 1.97 x 10"^ and eth = 1.56 x 10~\ The boundary condition of the field at 
Xi = 10~^ is = 1.2539010, which is larger than the the analytic value tfi{xi) = 1.09850009 that comes from Eq. (|76|l . The 

derivative <p'{xi) is chosen to be the same as the analytic value. The left panel depicts in the region < r/rc < 10, while 
the right panel depicts ip in the region < r/vc < 2 with log scales in the vertical and horizontal axes. The black curve (a) 
and the dotted curve (b) correspond to the numerically integrated solution and the analytic field profile (|76 p - (|78|l . respectively. 
The numerical solution outside the body recovers the analytic field profile (|78|) . 



given by the corresponding analytic value, we find that it is possible to reproduce the analytic thin-shell solution ([75]) 
outside the body even for <I>c ~ 0.1. The need for the choice of larger </) partially comes from the overestimation of the 
field around r = ri, as was explained above. Moreover, since the pressure is underestimated in our linear expansion 
of <i>c, we need to choose values of (f) larger than the corresponding analytic values in order to prevent the field from 
entering the region (j) < 4>A- 

In Fig. |4]wc plot an example of the numerical solution for <I>c = 0.2, n — 2, Q = 1, Arc/vc = 0.1 and niArc ~ 20.0, 
together with the corresponding analytic field profile. We have used the boundary condition 93 (cc^ = 10~^) = 1.2539010, 
which is larger than the analytic value (p{xi = 10~^) = 1.09850009 estimated by Eq. ([75)1 . We note again that the 
resulting field profile is sensitive to the change of boundary conditions. As can be seen from the right panel of Fig. [H 
the derivative 4>'{r) is negative in the region < r/vc < 0.69. The field grows for increasing r in the region r/rc > 0.69 
so that it enters the thin-shell regime for r/r^ > 0.9. The left panel of Fig. 2] shows that the numerical solution outside 
the body agrees well with the corresponding analytical solution. The solution asymptotically approaches the field 
value = 21.844. 

We have also carried out numerical simulations for other model parameters in the strong gravitational backgrounds. 
We find that thin-shell solutions are present for S 0-3, which marginally includes the case of neutron stars. When 
$c ^ 0.3, however, the field continues to evolve toward smaller (j) and overshoots the effective potential maximum at 
(j) = (pA (i-e., iy9 = 1) unless the boundary condition around the centre of the body is chosen to be ^ 1- The 

evolution of the field is typically followed by the rapid roll-down along the potential toward the singularity at = 
(as in the numerical simulations of Kobayashi and Maeda [2^ for the /(i?) dark energy model of Starobinsky [sot ). 
Since the ratio Pm (r) / pa is of the order of around the centre of the body, the pressure force is so strong that the 
field typically overshoots the effective potential maximum in such cases. We stress here that in strong gravitational 
backgrounds with $c = a separate analysis is required without recourse to the analytic solutions derived here 

which are valid only in the regimes with <i>c < 0(0.1). 

Finally we note that the distance ^3 at which (l)'{r^) = gets smaller for decreasing niA (see Table I). This may 
suggest that it is possible to avoid the overshooting of the field by choosing smaller values of rriA- However, the 
parameter Ar^/rc needs to satisfy the conditions Ar^jr^ <C 1 and Ardr^ 3> 1/toa?'c- This implies that we can not 
choose the values of to^Tc that are smaller than the order of 10. Thus when Arc/j'c ^ O.l and rriArc ^ 10, it is 
typically difficult to obtain thin-shell solutions for (f>c 0.3, whereas thin-shell solutions arc present for <f>c ^ 0.3. 



16 



VI. CONCLUSIONS 

111 this paper we have studied the behaviour of the chameleon scalar field (j) in the relativistic gravitational back- 
ground of the spherically symmetric space time. The gravitational potentials $ and \E' are found analytically under 
the conditions that the density of the central compact object is constant and that the energy density of the chameleon 
field is much smaller than that of the matter. Using the gravitational potential $c at the surface of the body as a 
linear expansion parameter we have derived the scalar-field equation (|17p . 

The solutions to the field equations can be obtained by considering the perturbation 54> about the corresponding 
solution 00 in the Minkowski background. In the region < r < ri the field exists around the minimum of the 
effective potential V^ff (0) = V{4)) + QpA<l> inside the body. The thin-shell case corresponds to settings in which ri is 
close to the radius rc of the body. In the region ri < r < rc the coupling term Qpa4> dominates over the effective 
potential, which leads to rapid changes in the field. Using linear expansions in terms of and (50 we have derived 
the solutions of the field equation in the regions < r < ri and ri < r < Tc, see Eqs. (PU)) and Outside the body, 
the kinetic energy of the field dominates over its potential energy, so that the approximate solution in this region is 
given by Eq. ((551) . 

We have matched the three solutions at the distances ri and subject to boundary conditions (pS)) and have 
derived the analytical thin-shell field profile given by Eqs. ([55)) - ([55)1 . In discussing the analytical field profile, we have 
considered the case Q > for simplicity. Compared to the result of the Minkowski spacetime, the field (j) around the 
centre of the body is shifted due to the presence of a relativistic pressure. For larger values of $c and ttiat'c, the field 
derivative (t>'{r) becomes negative in the region < r < ra (< ri). For values of r > rs, 0'(r) becomes positive and 
the field 0(r) begins to grow with increasing r. As long as the condition /S.rc/rc ^ l/(m^rc) is satisfied, the field 
acquires sufficient kinetic energy in the thin-shell regime in order to climb up the potential hill outside the body. 

For the class of potentials V{(t>) = Af''+"0~" we have carried out numerical simulations by using the information 
provided by the analytic field profile in order to set the boundary conditions around the centre of the body. In the 
Minkowski background (<I>c = 0) the thin-shell field profile outside the body can be recovered numerically by choosing 
the boundary condition of the field to be larger than the corresponding analytic value. The reason for this comes 
from the fact that the analytic solution overestimates the field value in the region r'l < r < by neglecting the term 
U^ relative to QpA- 

In the relativistic gravitational backgrounds with <f>c ^ 0.3 we have also confirmed the presence of thin-shell solutions 
numerically. While there exists a region in which (j)' {r) is negative inside the body, it is possible to realize thin-shell 
solutions if the derivative (j)' {r) changes sign at a distance r = smaller than ri. For larger $c the distance ^3 tends 
to increase so that the effect of the relativistic pressure is stronger inside the body. We note that our analysis does not 
cover the case of extremely stong gravitational backgrounds with of the order of unity. This requires a separate 
detailed analysis which incorporates the formation of black holes. 

Finally we note that realistic stars have densities pAix) that globally decrease as a function of r. It would be 
expected that this decreasing density may work as a counter term to the relativistic pressure around the centre of the 
body [see Eq. ([TT)) ]. It would be of interest to see whether thin-shell solutions arc present in such realistic cases with 
strong gravitational backgrounds. We shall return to this question in future. 
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